CERN-PH-TH/12-184, RM3-TH/12-12, DFPD-2012/TH/7, CP3-Origins-2012-018, DIAS-2012-19, TUM-HEP-843/12 



Repressing Anarchy in Neutrino Mass Textures 

Guido Altarelli a ^,* Ferruccio Feruglio c) , f Isabella Masina and Luca Merlo f^)% 

a ' Dipartimento di Fisica 'E. Amaldi', Universitd di Roma Tre, 
INFN Sezione di Roma Tre, 1-00146 Rome, Italy 
b ' CERN, Department of Physics, Theory Division, CH-1211 Geneva 23, Switzerland 
°> Dipartimento di Fisica e Astronomia 'G. Galilei', Universitd, di Padova, 
INFN Sezione di Padova, Via Marzolo 8, 1-35131 Padua, Italy 

d ' Dipartimento di Fisica dell'Universita di Ferrara, 
INFN Sezione di Ferrara, Via Saragat 1, 1-44100 Ferrara, Italy 
e) CP 3 -Origins and DIAS, SDU University, 
Campusvej 55, DK-5230 Odense M, Denmark 
Physik- Department, Technische Universitat Miinchen, 
James-Franck-Strasse, D-85748 Garching, Germany 
9 ' TUM Institute for Advanced Study, Technische Universitat Miinchen, 
Lichtenbergstrasse 2a, D-85748 Garching, Germany 



The recent results that #13 is relatively large, of the order of the previous upper bound, and the 
indications of a sizable deviation of #23 from the maximal value are in agreement with the predictions 
of Anarchy in the lepton sector. The quark and charged lepton hierarchies can then be reproduced 
in a SU(5) GUT context by attributing non- vanishing U(1)fn charges, different for each family, 
only to the SU(5) tenplet states. The fact that the observed mass hierarchies are stronger for up 
quarks than for down quarks and charged leptons supports this idea. As discussed in the past, in 
the flexible context of SU(5) <8> U(1)fn, different patterns of charges can be adopted going from 
Anarchy to various types of hierarchy. We revisit this approach by also considering new models and 
we compare all versions to the present data. As a result we confirm that, by relaxing the extreme 
ansatz of equal U(1)fn charges for all SU(5) pentaplets and singlets, better agreement with the 
data than for Anarchy is obtained without increasing the model complexity. We also present the 
distributions obtained in the different models for the Dirac CP-violating phase. Finally we discuss 
the relative merits of these simple models. 



I. INTRODUCTION 



Recently our knowledge of the neutrino mixing matrix has much improved with the rather precise 
measurement of #13 and the indication that O23 is not maximal [1, 2] (some constraints on the CP- 
violating phase in neutrino oscillations are also starting to emerge). The rather large measured value 
of #13 [3-6], close to the old CHOOZ bound [7] and to the Cabibbo angle, and the indication that 
#23 is not maximal both go in the direction of models based on Anarchy [8, 9], i.e. of no symmetry 
in the leptonic sector, only chance (this possibility has been recently reiterated, for example, in Ref. 
[10]). Anarchy can be formulated in a SU(5) (g> U(1)fn context by taking different Froggatt-Nielsen 
[11] charges only for the SU(5) tenplets (for example 10 ~ (a, 6,0), where a > b > is the charge 
of the first generation, b of the second, zero of the third) while no charge differences appear in the 
5 (e. g. 5 ~ (0,0,0)). If not explicitly stated, the Higgs fields are taken neutral under U(1)fn- 
The SU(5) generators act vertically inside one generation, whereas the U(1)fn charges are different 
horizontally from one generation to the other. If, for a given interaction vertex, the U(1)fn charges do 
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not add to zero, the vertex is forbidden in the symmetric limit. However, the U(1)fn symmetry (that 
we can assume to be a gauge symmetry 1 ) is spontaneously broken by the VEVs Vt of a number of 
flavon fields with non- vanishing charges and GUT-scale masses. Then a forbidden coupling is rescued, 
but is suppressed by powers of the small parameters A = Vf/M, with M a large mass, with the 
exponents larger for larger charge mismatch. Thus the charges fix the powers of A, hence the degree of 
suppression of all elements of mass matrices, while arbitrary coefficients c^- of order 1 in each entry of 
mass matrices are left unspecified (so that the number of parameters exceeds the number of observable 
quantities). A random selection of these Cy parameters leads to distributions of resulting values for the 
measurable quantities. For Anarchy (A) the mass matrices in the leptonic sector are totally random; 
on the contrary, in the presence of non-vanishing charges different entries carry different powers of the 
order parameter and thus some hierarchies are enforced. There are many variants of these models: 
fermion charges can all be non-negative with only negatively charged flavons, or there can be fermion 
charges of different signs with either flavons of both charges or only flavons of one charge. In models 
with no SeeSaw, the 5 charges completely fix the hierarchies (or Anarchy, if the case) in the neutrino 
mass matrix. If Right-Handed (RH) neutrinos are added, they transform as SU (5) singlets and can in 
principle carry U(1)fn charges, which also must be all equal in the Anarchy case. With RH neutrinos 
the SeeSaw mechanism can take place and the resulting phenomenology is modified. 

Anarchy and its variants, all sharing the dominance of ramdomness in the lepton sector, are to be 
confronted with models based on discrete flavour groups. These more ambitious models are motivated 
by the fact that the data suggest some special mixing patterns as good first approximations like Tri- 
Bimaximal (TB) or Golden Ratio (GR) or Bi-Maximal (BM) mixing, for example. The corresponding 
mixing matrices all have sin 2 823 = 1/2, sin 2 #13 = 0, values that are good approximations to the 
data (although less so since the most recent data), and differ by the value of the solar angle sin 2 812- 
The observed sin 2 #12, the best measured mixing angle, is very close, from below, to the so called 
Tri-Bimaximal (TB) value [12-16] of sin 2 8 12 = 1/3. Alternatively, it is also very close, from above, 
to the Golden Ratio (GR) value [17-20] sin 2 8 12 = l/V5</> ~ 0.276, where <f> = (1 + V$)/2 is the GR 
(for a different connection to the GR, see Refs. [21, 22]). On a different perspective, one has also 
considered models with Bi-Maximal (BM) mixing [23-26], where sin 2 812 — 1/2, i.e. also maximal, 
as the neutrino mixing matrix before diagonalization of charged leptons. One can think of models 
where a suitable symmetry enforces BM mixing in the neutrino sector at leading order (LO) and the 
necessary, rather large, corrective terms to 9±2 arise from the diagonalization of the charged lepton 
mass matrices [27-54]. Thus, if one or the other of these coincidences is taken seriously, models where 
TB or GR or BM mixing is naturally predicted provide a good first approximation (but these hints 
cannot all be relevant and it is well possible that none is). The corresponding mixing matrices have 
the form of rotations with fixed special angles. Thus one is naturally led to discrete flavour groups. 
Models based on discrete flavour symmetries, like A4 or S4, have been proposed in this context and 
widely studied [55-61]. In these models the starting Leading Order (LO) approximation is completely 
fixed (no chance), but the Next to LO (NLO) corrections still introduce a number of undetermined 
parameters, although in general much less numerous than for U(1)fn models. These models are 
therefore more predictive and typically, in each model, one obtains relations among the departures of 



1 Gauge anomalies can be cancelled adding a set of heavy additional fermions, vector-like under SU(5) and chiral under 
U(1)fn- They have no impact on the present discussion. 
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the three mixing angles from the LO patterns, restrictions on the CP violation phase S, mass sum rules 
among the neutrino mass eigenvalues, definite ranges for the neutrinoless-double-beta decay effective 
Majorana mass and so on. 

The aim of this note on U(1)fn models is to update, on the basis of the present more precise data, 
our old analysis [62] that shows that, even if one accepts a mainly chaotic approach to lepton mixing, 
the Anarchy ansatz is too extreme and that suitable differences of U(1)fn charges, which must in 
any case be present for tenplets, if also introduced within pentaplets and singlets, actually lead to 
distributions that are in much better agreement with the data with the same number of random 
parameters. In fact Anarchy can be improved by implementing mechanisms that enforce the relative 
smallness of $13 and of r = Am^ olar / Anr^ tm . The first goal can be achieved by restricting Anarchy 
only to the fj, — r (or 2-3) sector, by taking the pentaplet charges as 5 ~ (c, 0,0), c > 0, ("/i — r 
Anarchy", A MT ) [63, 64]. Both improvements can be realised by taking 5 ~ (c, d, 0), c > d > 0, (this 
"hierarchical", H, pattern was not considered in Rcf. [62]). In each case the tenplet charges can 
be readjusted in order to maintain the correct ratios of charged fermion masses. We only consider 
models with normal hierarchy (NH) here, because, as shown in Ref. [62], in this framework models 
with inverse hierarchy (IH) tend to favour a solar angle close to maximal. In models with no SeeSaw, 
the 5 charges completely fix the hierarchies (or the Anarchy) in the neutrino mass matrix, through the 
dimension-5 Weinberg operator, m v = ^g^g H u H u /M. The distributions arising from the models 
A, A^t and H with no SeeSaw can directly be obtained and compared with the data and we shall 
see to which extent H is better than A^ T which, in turn, is better than A. If RH neutrinos are 
added, they transform as SU(5) singlets and can in principle carry U(1)fn charges, which also are all 
equal in the Anarchy case. With RH neutrinos the SeeSaw mechanism can take place and the resulting 
phenomenology is modified. It is easy to show that models with all non-negative charges and one single 
flavon have particularly simple factorization properties [65]. In particular, in the SeeSaw expression 
for the light neutrino mass matrix m v — mjj M~ l mc, with mj and M denoting the neutrino Dirac 
and Majorana mass matrices, respectively, the dependence on the RH charges drops out in this case 
and only that from the 5 remains. In these simplest models the only difference between the version 
with and without SeeSaw, for each model A, A^ T and H is that the extraction procedure for the 
random numbers is different: in the no SeeSaw version the entries of the neutrino mass matrix m„ 
are directly generated while in the SeeSaw case the extraction is done for mr> and M and then m„ 
is derived by the SeeSaw formula. For example, in Anarchy models the smallness of r is to some 
extent reproduced by the spreading of the mass distribution resulting from the product of three mass 
matrices. Models with naturally large 23 neutrino mass splittings (so that r is small) are obtained 
if we allow negative charges and, at the same time, either introduce flavons of opposite charges or 
allow that matrix elements with overall negative charge are vanishing. For example, one can take 
5 ~ (c, 0, 0) like in "/i — r" Anarchy and 1 ~ (e, — e, 0), e > 0, with two flavons of opposite charges, 
with equal VEV, and SeeSaw (we denote this model as Pseudo /zr- Anarchy, "PA^t"). The "lopsided" 
structure of 5 ~ (c, 0, 0) results in naturally small 23 subdeterminant in the neutrino mass matrix 
after SeeSaw and to r naturally small. 

In the following we discuss in more detail the models, the procedure of extraction of the random 
coefficients and the value of the expansion parameter A that maximizes the success rate for each model. 
We then discuss the mass distributions and mixing angles that we obtain and, finally, we compare 
these distributions with the data. The conclusion is that the most effective model is H in the no 
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SeeSaw case and PA^t in the SeeSaw case, while A^ T and A follow in the order and are much less 
successful. 

II. MODELS AND RESULTS 

In the following analysis we adopt the results of the fit of Ref. [1] (see also Ref. [2]). The 2(3)er 
ranges in the case of NH are: 

7.15(6.99) x KT 5 eV 2 < Am 2 solar < 8.00(8.18) x 10" 5 eV 2 
2.26(2.15) x 10~ 3 eV 2 < Am 2 atm < 2.58(2.66) x 10~ 3 eV 2 

0.0181(0.0149) < sin 2 6 13 < 0.0311(0.0344) (1) 

0.275(0.259) < sin 2 6» 12 < 0.342(0.359) 

0.350(0.330) < sin 2 6*23 < 0.475(0.638) 

where Am 2 o;ar = m\ — m\ and Am 2 tm = m§ — (m 2 + m\)/2. 

We consider models with different patterns - Anarchy, genuine or Pscudo fir- Anarchy, and hierarchy 
- induced by a U(l) flavour symmetry [11]: we present the transformation properties of all the fields 
in table I, in a notation that corresponds to the SU (5) GUT embedding. In the non-SeeSaw case, the 
neutrino mass matrix m v is generated via the effective Weinberg operator \E r g^\E , 5 H u H U JM, In the 
SeeSaw case the flavour charges determine the Dirac and RH Majorana mass matrices, mo and M, 
which give the effective neutrino mass m v = rn^M^m^, at low energy. If the RH neutrino charges 
all have the same sign and there is a single flavon, it is known that the structure of m„ in powers of 
A is the same as for the non-SeeSaw case [66]. The coefficients in front of A are randomly generated 
complex numbers c = | c| e*^ . In the spirit of the U(l) flavour symmetry, |c| = 0(1) while the phase 
4> c is arbitrary. 

We report in Fig. 1 the per cent success probability for each model as a function of the expansion 
parameter A. We defined P — n k/n to t, where n to t is the total number of randomly generated models 
(typically larger than 10 6 ) and n a k is the number of models consistent with the 2a ranges of Eq. (1) 
[1]. The relative error on P is estimated to be 0(l/-i/n o fe). The thickness of the lines indicates 
this statistical error. Since the success rate depends on the window selected to extract the random 
coefficients, to estimate the ambiguity from this effect we flatly generated |c| in the interval [0.5, 2] 
(solid shaded) and also in the interval [0.8, 1.2] (dashed shaded). The phases <fi c have been chosen to 
be flatly distributed in [0, 2ir]. 

The flavour charges of in table I have been chosen in order to reproduce the mass hierarchies of 
the charged fermions for the value of A that for each model maximizes the success rates. We scanned 
for V^io integer charges (a, 6, 0) with < b < a < 10 and b < 5. For each choice of charges we compared 
the distributions of the six mass ratios m e /m^, m^/rUr, md/m s , m s /m(,, m u /m c and m c /m t and of 
the three CKM matrix elements V us , V u b and V c b with the corresponding 3tr experimental interval, 
rcnormalized at the GUT scale assuming tan/3 = 10 (see for instance Ref. [67]). We extracted the 
most successful charges by asking that the experimental interval of each parameter overlaps with the 
la region of the theoretical distribution. Only for the PA^ T model, this requirement is fulfilled for all 
the nine observables. For all the other models, there is a tension between the choices which fit m e /m„, 
which favour large values of a — b, and the choices that reproduce V us , which require a — b small. In 
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Model 


#10 


#5 




Anarchy (^4) 


(3,2,0) 


(0,0,0) 


(0,0,0) 


/it- Anarchy (j4 mt ) 


(3,2,0) 


(1,0,0) 


(2,1,0) 


Pseudo fir- Anarchy (PA^t) 


(5,3,0) 


(2,0,0) 


(1,-1,0) 


Hierarchy (H) 


(5,3,0) 


(2,1,0) 


(2,1,0) 



TABLE I: Models and their flavour charges suitable for an implementation in a Supersymmetric (SUSY) 
SU(5) Cg> C/(1)fn GUT. The flavon charge is —1. The charges of the #io have been chosen to reproduce the 
mass hierarchies of the charged fermions, for the values of A that maximize the success rates for each model 
(see text). The Higgs H Ui d charges have all been taken as vanishing in these models. 



these models, a tuning of the unknown order one parameters is needed in order to reproduce both 
these quantities. We gave our preference to the solutions which correctly reproduce V us . For a given 
model, the solution we found is not unique and other choices of the charges are equally successful, 
provided we change the size of the expansion parameter A: indeed, redefining A — > A 1 / 9 , a — > qa 
and b — ?• q b, where q is a positive number, the structure of the mass matrices that we will show in 
the following does not change. This observation alleviates the otherwise arbitrary choice of integer 
charges. In table I we list a representative set of charges that passed our test. Choosing a vanishing 
charge for both the 10 and the 5 SU(5) representations of the third generation have no impact on our 
selection procedure, that relies on quantities that are only sensitive to charge differences. The nominal 
values of the charges listed in table I require a large value of tan /3, since the top and bottom Yukawa 
couplings are both expected to be of order one. Smaller values of tan/3 can be easily accommodated 
by allowing for a positive charge of H^. 

In Fig. 2 we display the probability distributions for r, sin# 13 , tan 2 # 12 and tan 2 # 2 3, fixing A at a 
representative value and selecting |c| S [0.5,2]. These distributions have been normalized so that the 
integrated probability is equal to unity. Notice that we use a logarithmic scale. The (green) shaded 
vertical regions refer to the experimental data at 2a [1]. 

We now discuss each model in turn. 

• Anarchy (A) [8-10]. Neglecting the randomly generated coefficients, the texture of the mass 
matrices for charged leptons and neutrinos (with or without SeeSaw), expressed in powers of A, 
reads explicitly 



/a 3 A 3 A 3 \ (l 1 1 N 

A 2 A 2 A 2 , ro„ = 1 1 1 | . (2) 

V 1 1 V V 1 h 



G 



0.010 



o 
o 



0.005 



Non-SeeSaw 




_ 




''/ H /\\ 























0.000 

0.1 0.2 0.3 0.4 



0.010 



0.5 



0.6 




FIG. 1: Per cent probability of success to fulfill the 2a ranges of Eq. (1) [1] as a function of A, without 
(left) and with SeeSaw (right). Solid (dashed) lines are obtained extracting |c| with a flat distribution in the 
interval [0.5, 2] ([0.8, 1.2]). Phases <f> c are flatly distributed in [0, 2n]. The thickness of the curves represent the 
statistical error, estimated as discussed in the text. In the left plot, the success rate for A has been multiplied 
by a factor of 10. 



The success rate for neutrino masses and mixing angles is independent of A. We then choose 
A = 0.2 — 0.3, which ensure a reasonable hierarchy for charged fermions according to the charges 
selected for the \l/io representation. As can be seen from Fig. 1, the success rate is quite small 
in both the no SeeSaw and SeeSaw cases: for the no SeeSaw case the success rate has been 
multiplied by a factor of 10 to facilitate its comparison with the other models. The reasons for 
such a modest performance can be understood by inspecting Fig. 2: the most severe problem is 
the prediction of a too large #13 and, in the no SeeSaw case, also of a too large value of r. In the 
SeeSaw case the latter problem is cured by the spreading of neutrino mass eigenvalues produced 
by the product of three random matrix factors. As for the mixing angles fty the distributions 
are all similar and, with a logarithmic scale, appear peaked near 7r/4. 

fir- Anarchy (A MT ) [63, 64, 68]. In this case only the /it block of m„ is anarchical 

/a 4 a 3 a 3 \ /a 2 a a\ 

m t = A 3 A 2 A 2 , m u = A 1 1 (3) 

\X ll) \X 1 l) 

and the success rate is maximized for A = 0.2 and A = 0.25 for the no SeeSaw and SeeSaw cases 
respectively. In both cases the performance of Ap T is better than A. The main problem of the 
A^t model is the prediction of a too small 612 and, in the no SeeSaw case, also a too large value 
for r. If by accident the 22 matrix element of m v is numerically of order A then # 12 ~ 0(1) and 
i/r ~ O(A): with a single fine tuning one fixes both problems. 

Pseudo [it- Anarchy (PA^ T ) [62]. This is a SeeSaw model with two flavons of opposite charges 
and equal VEVs. The 2 and 3 entries of the pentaplets have the same charges, but the 1 and 
2 RH neutrinos have opposite charges. The result is that m„ displays an apparently anarchical 
23 sector, 




m t = A 5 A 3 A 3 , m D = A A A , M = 1 A 2 A , m v = A 2 1 1 , (4) 
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but the associated coefficients automatically induce a suppression of the /it determinant 2 , which 
is desirable to justify the smallness of r while <9 12 ~ The success rate for this model is 

maximized for A = 0.35 — 0.4. For such values the distributions of r and tan 2 6*12 in Fig. 2 are 
indeed nearly centered in the experimental range. The mixing angle #23 is instead naturally 
maximal and its distribution is indistinguishable with respect to A and models. Notice that 
the PA^ model emerged as favorite in the 2005 update of the analysis of Ref. [62] . 

• Hierarchy (H). Both for the SeeSaw and no SeeSaw cases, the charged lepton and neutrino mass 
matrices read 



The success rate is maximized for A = 0.4 and A = 0.45 for the non-SeeSaw and SeeSaw cases, 
respectively. For non-SeeSaw one finds a successful model every 10, 000 trials, that is a factor of 
100 better than the Anarchy texture. For the SeeSaw the success rate is slightly lower, but still 
a factor of 10 better than Anarchy. This is mainly due to the hierarchical structure in both 12 
and 23 sectors, that ensures a small r ~ A 4 and tan 2 f? 12 ~ tan 2 #23 ~ sin #13 ~ A 2 . As can be 
seen from Fig. 2, with A = 0.4, the maxima of the distributions of these observables are nicely 
close to their experimentally allowed range at 2a. Notice that the distributions of 6*12 and 6*23 
are similar. In particular, #23 is peaked at a slightly smaller value than its present experimental 
2cr range. Despite this, the rate of success of 6*23 equals those of the three previously studied 
versions of anarchical models, namely A, A^ T and PA^. 

We can now directly compare the four models. As from Fig. 1, H is the best performing model 
for the non-SeeSaw case, for values of A larger than about 0.3, while for smaller values, A^ T has the 
best success rate. This suggests that a moderate hierarchy could likely be realized in the neutrino 
sector. For the SeeSaw case, the performances of H and PA^ T are almost equivalent, although PA^r 
is slightly better. 



It is also interesting to study the distribution of the Dirac CP-violating phase 5 in the models 
considered. At present, there is just a very mild la preference for 8 ~ it in the NH case [1]. 
In order to extract the phase (5, we consider the following combination 

that is an invariant under phase transformation of the fields. This is only one of the possible invariants 
that can be considered (notice that the imaginary part of Z7 e 3 U* 2 U* 3 corresponds to the Jarlskog 





III. DIRAC CP PHASE 



2 Without the RH neutrinos, PA^ T model corresponds to the A^ T , as indeed this suppression mechanism for the 
determinant does not hold. 
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FIG. 2: Probability distributions of r, sin #13, tan 2 #12, tan 2 #23 without (left column) and with (right column) 
SeeSaw. The modulus (argument) of the complex random coefficients has been generated in the interval 
[0.5, 2] ([0, 2n]) with a flat distribution. For A and A^ lT we considered A — 0.2, for H and PA^ we considered 
A = 0.4. The shaded vertical band emphasizes the experimental 2a window according to [1]. 
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FIG. 3: Probability distributions of S without (left column) and with (right column) SeeSaw. The modulus 
(argument) of the complex random coefficients has been generated in the interval [0.5,2] ([0, 2tt]) with a flat 
distribution. For A and A^ r we considered A = 0.2, for H and PA^ we considered A = 0.4. Note that we 
considered S 6 [— n,n]. 

invariant [69]). Adopting the usual PDG parameterisation of the PMNS matrix, we get 



T If/, 



,.,2 



COS 023 COS 012 e 



-iS 



sin 6*23 sint/i3 Sin #12 



(7) 



The distributions of 5 are shown in Fig. 3. As one expects, for the Anarchy model the distribution 
of 5 is completely flat. On the other hand, for H and A^ T there is a mild preference for a vanishing 
value of 5, while for PA^ T this preference is even weaker. 



IV. CONCLUSION 



Over the years there has been a continuous progress in the measurement of neutrino mixing angles 
culminating recently with the determination of a relatively large value of 6*13 and with the indication 
that 023 is not maximal (some hints that cos 8 < 0, with S being the Dirac CP-violating phase in 
neutrino oscillations, are also emerging). In spite of this remarkable experimental progress the data 
can still be reproduced by a wide range of theoretical models. At one extreme we have models where the 
assumed dynamics is minimal and the dominant ingredient is pure chance (Anarchy and its variants) 
and, at the other extreme, models with a high level of underlying symmetry, like, for example, those 
based on discrete non-Abelian symmetries (which start at LO with TB or BM mixing). The large 
value of 0i3 and the departure of #23 from maximal both go in the direction of Anarchy and move away 
from the TB or BM limits, where $13 = and 023 is maximal. In this note we have made a reappraisal 
of Anarchy, given the new experimental results. To make connection with quark masses and mixing 
we have adopted the (SUSY) SU(5) <g> ?7(1)fn GUT framework. The Anarchy prototype model has 
only tenplet charge differences (among the 3 generations) that are non-vanishing, while all pentaplet 
and singlet charge differences are taken as vanishing. Here we confirm on the most recent data that 
the Anarchy ansatz is too extreme. In fact, suitable differences of U(1)fn charges, if also introduced 
within pentaplets and singlets, lead to distributions that are in much better agreement with the data 
with the same number of random parameters as for Anarchy. The hierarchy of quark masses and 
mixing and of charged lepton masses in all cases impose a hierarchy defining parameter of the order 
of Ac = sin 8c, with 9c being the Cabibbo angle. The weak points of Anarchy (A) are that with this 
ansatz all mixing angles should be of the same order, so that the relative smallness of 6*13 ~ O(Xc) 
is not automatic. Similarly the smallness of r = Am^ olar / Am^ tm is not easily reproduced: with no 
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SeeSaw r is of 0(1), while in the SeeSaw version of Anarchy the problem is only partially alleviated by 
the spreading of the neutrino mass distributions that follows from the product of three matrix factors 
in the SeeSaw formula. An advantage is already obtained if Anarchy is only restricted to the 23 sector 
of leptons as in the model. In this case, with or without SeeSaw, #13 is naturally suppressed and, 
with a single fine tuning one gets both #12 large and r small. Actually we have shown that, in the 
no SeeSaw case, a very good performance is observed in the new H model, where the Anarchy is also 
abandoned in the 23 sector. In the H model, by taking a relatively large order parameter, one can 
reproduce the correct size for all mixing angles and mass ratios. In the SeeSaw case, we have shown 
that the freedom of adopting RH neutrino charges of both signs, as in the PA^ r model, can be used 
to obtain a completely natural model where all small quantities are suppressed by the appropriate 
power of A. In this model a lopsided Dirac mass matrix is combined with a generic Majorana matrix 
to produce a neutrino mass matrix where the 23 subdeterminant is suppressed and thus r is naturally 
small with unsuppressed 6*23- In addition also #i 2 is large, while 6>i 3 is suppressed. We stress again 
that the number of random parameters is the same in all these models: one coefficient of 0(1) for 
every matrix element. Moreover, with an appropriate choice of charges, it is not only possible to 
reproduce the charged fermion hierarchies and the quark mixing, but also the order of magnitude of 
all small observed parameters can be naturally guaranteed. Thus finally we agree that models based 
on chance arc still perfectly viable but wc consider Anarchy a too drastic and unjustified abstraction. 
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